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Abstract 

A recent CMS analysis has reported the observation of an excess in the invariant 
mass distribution of the opposite-sign same-flavour lepton pair, which can be inter¬ 
preted as a kinematic edge due to new physics. Using collider simulation tools, we 
recast relevant LHC search results reported by ATLAS and CMS collaborations in 
order to determine constraints on supersymmetric models that could produce the 
observed features. In particular, we focus on models involving cascade decays of 
light-flavour squarks and sbottoms. We find no favourable supersymmetry scenario 
within our exploration that could explain the origin of the excess when other LHC 
constraints are taken into account. 


^ philipp. grothaus@kcl .ac.uk 
^ liew@hep-th. phys .s .u-tokyo .ac.jp 
^kazuki.sakurai@kcl.ac.uk 



1 Introduction 


The search for supersymmetry (SUSY) as an extension of the Standard Model (SM) is 
one major target of the Large Hadron Collider (LHC) physics program. However, we have 
thus far found no definitive evidence of SUSY based on the first run of the LHC despite 
dedicated searches on many fronts. Even so, there is an analysis presented by the CMS 
collaboration that could be showing the first signs of SUSY BE]- In the analysis, two 
leptons, jets and missing energy are looked for in the final states. It is found that there is 
an excess (130^4g events in the “central” region) in the invariant mass distribution of the 
opposite-sign same-flavour (OSSF) lepton pair, corresponding to a signihcance of 2.6 a. 

The excess of the signal is htted kinematically as a triangular-shape edge at mu = 
78.7 ± 1.4 GeV. Such a kinematic edge is a characteristic signal of SUSY, where a SUSY 
particle undergoes a two-stage two-body decay. The kinematic edge formed by a pair of 
leptons can be interpreted as the cascade decay of a neutralino: y® —>■ 

(on-shell slepton decay) [3] . It is also possible to interpret the edge as a three-body decay 
signal of a neutralino, y® where the lepton pair is produced via an off-shell Z 

(off-shell Z decay). The shape of the edge would be more rounded as compared to the two- 
stage two-body decay, but as shown in the original CMS analysis, the three-body decay 
still provides a good fit. The direct production of X 2 too small to reproduce the dilepton 
excess, however its production can be boosted if coloured sparticles subsequently decay 
into X 2 - The explanation of the dilepton excess in terms of coloured sparticles is consistent 
with the CMS analysis, since events with jet multiplicity are selected and counted. 

In this work we perform a detailed study on the possibility of explaining the dilep¬ 
ton excess with several SUSY models taking into account a comprehensive list of LHC 
constraints from a number of ATLAS and CMS direct SUSY searches. In order to accu¬ 
rately estimate the LHC constraints and simulate many analyses systematically, we use 
the automated simulation tool Atom |1]. We take a bottom-up approach by considering 
simplihed SUSY models with minimal content of particles at low energy to reproduce the 
excess optimally. As will be discussed in the following sections, light-flavour squarks and 
sbottoms are potential candidates that satisfy these criteria. Some of these models have 
already been studied in earlier works Here, we will confront our simplihed models 

with various direct SUSY search constraints such that their viability is tested in great 
detail. We will show that the light-havour squarks and sbottom models we consider in 
this paper are strongly constrained when providing a large enough contribution to the 
dilepton excess. 

Our paper is organised as follows. In the next section, we describe the selection criteria 
of the CMS dilepton analysis. In section we consider SUSY models that can possibly 
reproduce the required features of the dilepton edge. In section we describe the proce¬ 
dure of our simulation and analysis. We discuss our results and their interpretations in 
Section m Conclusions are drawn in Section |6l 

^See [7] for a non-SUSY interpretation of the observed excess. 
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2 CMS dilepton analysis 

CMS reported an excess of events in the dilepton plus missing energy channel DEI in the 
8 TeV, 19.4 fb“^ data. The analysis requires an OSSF lepton pair with pt > 20 GeV. It 
also requires > 2 jets with pt > 40 GeV and > 150 GeV or > 3 jets with pt > 40 
GeV and > 100 GeV. The excess is observed in the central region where both leptons 
satisfy l^^iepl < 1.4. It exhibits an edge in the dilepton invariant mass distribution around 
= 78 GeV. The counting experiment in the G [20, 70] GeV region shows an excess 
of ~ 130 events over the Standard Model expectation, which corresponds to a standard 
deviation of 2.6 a. 


3 SUSY interpretations of the dilepton edge 

In this paper we consider simplihed SUSY models that capture the essence needed for 
explaining the observed dilepton excess. Generalizations of SUSY models given here are 
straightforward. 

It is known that the OSSF dilepton pair in the decay of the second lightest neutralino 
X 2 via on-shell slepton and off-shell Z exhibit an edge-like shape at 


^edge 



^edge = : 


X 2 

xl ^ 


( 1 ) 

( 2 ) 


respectively, in the mgg distribution. In order to obtain a large enough production cross 
section to £t the excess and to have > 2 high px jets required in the event selection, we 
consider production of coloured SUSY particles, which may subsequently decay into X 2 - 
In this paper we consider two scenarios: light-flavour squark and sbottom production. 


3.1 Squark scenarios 

In the squark scenario, we consider the production of pairs of light-flavour squarks. This 
scenario assumes the hrst two generations of squarks (both left and right-handed) to be 
mass degenerate and within the reach of the LHG, whilst the third generation squarks 
and gluinos are decoupled. We also assume that the second lightest neutralino is mostly 
Wino-like or an admixture of Wino and Higgsinos and the lightest neutralino is mostly 
Bino-like. In this setup the lighter chargino, Xi: is naturally introduced as a SU{2)l 
partner of the X 2 their masses have to be close to each other. Since the right-handed 
squarks do not couple to the Wino and only very weakly couple to the Higgsinos, they 
decay predominantly into a quark and the Xi^ whereas the left-handed squarks can decay 
to either x^, X 2 or Xi depending on the Wino-Higgsino mixing in the Xi" ^rid y®- 
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We consider two models according to the xf and X2 decay modes. The hrst model is 
the on-shell slepton model, where we assnme the right-handed selection and smnon in the 
low energy spectrnm so that the X 2 decays predominantly into an OSSF lepton pair and 
the via the on-shell slepton. We deconple the left-handed slepton donblets, to 

maximise the signal rate, otherwise the X 2 conld also decay into a pair of nentrinos and 
the Xi via the on-shell 

Any models that lead to mnlti-lepton hnal states are severely constrained by the mnlti- 
lepton pins missing energy searches HD]. To avoid these constraints a large branching ratio 
of the qi —)■ qxi mode is necessary in this model. We assnme 


BR{qL ^ q + xf/xVxi) = 10/5/85 %. (3) 

This can be achieved if xf and X 2 have large Higgsino components becanse the sqnarks 
conple to the Higgsinos with small Ynkawa conplings. We will see in section 5.1 that onr 
conclnsion is not sensitive to variations of the branching ratios. 

In the sqnark with on-shell slepton model we then have the following cascade decays 


to 0 

H- 


5 %, 

Ql —t qxf —t qui^ — )■ qi^i^Xi ■ 

10 %, 


qL qxi ■ 

85%, 


qR qx[ ■ 

100% 


If one of the pair prodnced sqnarks nndergoes the hrst decay chain, the hnal state may 
contain an OSSF dilepton pins two energetic jets, and snch events can contribnte to the 
CMS dilepton excess. 

The second model we consider in this paper is the off-shell Z model, where the X 2 
decays via the oh-shell Z into an OSSF dilepton pair and the Unlike in the on-shell 
slepton model we here need a large branching ratio of q^ —t qx 2 snch that the small leptonic 
branching ratio of the oh-shell Z into electrons and mnons (abont 6 %) is compensated. 
We assnme 


BR{qL ^ q + xf/xl/xl) = 66/33/1 %, 


( 4 ) 


which can be realised by assnming xf and X 2 are Wino-like. In the sqnark with oh-shell 
Z model we have the following decay chains 


h qxl qffXi (via z*) 
Ql qxf qffxi (via W*) 
Ql qxl 
qr qxi 


33%, 
66 %, 
1 %, 
100 %. 


X2- 


®In our setup the X 2 decays predominantly into muon pairs through the Higgsino component of the 
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BR{qn j + Xi) = 100 % BR{qR j + Xi) = 100 % 

BR{qL ^ j + x^/xVxl) = 10/5/85 % BRiq^ ^ j + xf/x?/x?) = 66/33/1 % 



Figure 1: Decay chains of squark scenarios {left: on-shell slepton model, right: off-shell Z 
model). 


The signal events can be obtained if one of the pair produced squarks undergoes the hrst 
decay chain and the X 2 decays via the Z* into the dilepton pair and the x^. A schematic 
picture of the squark scenarios is shown in Fig. [Tj 

3.2 Sbottom scenarios 

Another way of interpreting the CMS dilepton excess is to assume that the observed 
dileptons in the excessive events come from cascade decays of bottom squarks. Unlike in 
the squark scenario, the decay mode to charginos, bi —)■ tx^, is kinematically forbidden 
if < rrit + m-±. We consider this case because the decay mode to charginos is more 
constrained due to emergence of top quarks. Similarly to the squark scenario we consider 
on-shell slepton and off-shell Z models according to the y® decay mode. 

In the on-shell slepton scenario the X 2 may decay either via a right-handed charged 
slepton or a left-handed charged slepton and sneutrino. We will treat these two cases 
separately in our analysis. 

If the mediating slepton is right-handed, X 2 predominantly decays into two charged 
leptons and and the events tend to have more than two leptons in the hnal state. Such 
models are severely constrained by the multi-lepton plus missing energy searches as we 
have previously discussed. To avoid these constraints, we assume 70 % of sbottoms decay 
into a bottom quark and a y? and the rest of sbottoms decay into a bottom quark and a 
y®- This situation can be achieved if y® is Wino-like and 61 has a large component of bn- 
We have the following decay chains in the sbottom with on-shell slepton model. 

bi bxl bi^i^xl : 30 %, 

bi ^ : 70 %. 
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Bii(6i^6 + X^/x?) = 50/50% 


+ = 30 / 70 % 


bi . h 



Figure 2: Decay chains of on-shell slepton mediated sbottom scenarios {left: left-handed 
slepton model, right: right-handed slepton model). 


In the case where the mediating slepton is left-handed, sneutrinos are introduced as 
SU(2) partners of charged sleptons. We assume sneutrinos and charged sleptons are mass 
degenerate and X 2 decays democratically into charged sleptons and sneutrinos. 


h bxl bi^i^ bi^i^Xi ■ 

25%, 

bi —)■ bx2 bui) —)■ buvx\ : 

25%, 

h bx\ : 

50%. 


The schematic picture of these cases is shown in Fig. 

Let us discuss the off-shell Z model for the sbottom scenario. Analogous to the squark 
with off-shell Z model, we need sbottoms to have a sizeable decay branching ratio to 
X 2 in order to have large enough dilepton event rates. One way to realise this situation 
is to have a Higgsino-like ^ mostly right-handed bi and to assume a large sbottom- 
bottom-Higgsino coupling due to a large tan//. It is shown in [6] that for tan// = 50, 
mi_^ ~ 330 GeV and a Higgsino mass parameter /r ~ 290 around 44% of sbottoms decay 
to the roughly mass-degenerate X 2 nnd Xa- This model point predicts about 1-a less 
events than the central £t without being excluded. In order to explore in more detail the 
parameter region that could possibly contribute to the excess, we expand the study of this 
scenario by varying the parameters Mi, /i and while hxing tan// = 50. The mass 
spectrum and particle decay branching ratio of this simplihed model are calculated using 
SPheno mu- 

Alternatively one can obtain a large branching ratio of sbottom decaying to X 2 by 
assuming bi is left-handed and X 2 Wino-like. Due to SU{2) gauge invariance a left- 
handed top squark, ti, is necessarily included in the low energy spectrum. For simplicity, 
we assume = mi^. We consider the following decay chains for the left-handed sbottom 
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Figure 3: Decay chains of off-shell Z mediated sbottom scenarios {left: right-handed 
sbottom model, right: left-handed sbottom model). 


with off-shell Z model. 


h -t Qxl QffXi (via Z*) : 100 %, 

ii Qxf Qff'Xi (via W*) : 100 %. 

A schematic picture of the off-shell Z sbottom scenarios is shown in Fig. 

The squark and sbottom scenarios are the priorities of this work, but let us also touch 
on the possibilities of explaining the dilepton excess with the remaining coloured sparticles 
in SUSY, namely gluino and stop. Gluinos can decay into X 2 via an intermediate squark, 
not much different from the squark or sbottom scenario other than a larger jet multiplicity. 
For stop, its decay into a top quark would lead to an extra lepton that plays no role at 
explaining the dilepton excess. It is not clear how gluino or stop could explain the dilepton 
excess without inducing additional jet or leptonic constraints, and hence we are not going 
to study these scenarios further in this work. 

4 The simulation setup 

In this section we describe our procedure to calculate the contribution to the CMS dilepton 
excess and the constraints from other ATLAS and CMS SUSY searches. 

The production cross section, cxprod, for light-flavour squarks is calculated using Prospino 
2 [11] with the gluino mass set to 3.5 TeV. For the sbottom cross section we use results 
from the LHC SUSY Cross Section Working Group based on [T2|. We create SLHA files 
of our simplihed models for the event generation and pass them to Pythia 6.4 [T3| to 
generate a total number of 10 ■ Uprod ■ with maximal 5 ■ 10®, events, where C = 19.4 fb“^ 
is the integrated luminosity at the CMS dilepton analysis. We then run Atom [1] on the 
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channel 

search for 

arXiv or CONF-ID 

refs 

2—6j -|- Of' -|- 


ATLAS-CONF-2013-047 

ra 



1405.7875 

ng 

2b + 0i + $T 

i, b 

1308.2631 

izni 

4j + li + $]T 

i 

ATLAS-CONF-2013-037 

IZB 

> 2j+ >li + $T 

q,g (1 or 2i) 

ATLAS-CONF-2013-062 

m 

2j -j- 2£ -j- ]^rp 

dilepton edge 

CMS-PAS-SUS-12-019 

BEI 

2j + + $T 

q,g,i,b (SS lepton) 

ATLAS-CONF-2013-007 

ra 

2j 2i f/rp 

i{2i) 

ATLAS-CONF-2013-048 

m 



1403.4853 

Ea 

2,3i fjj' 


1404.2500 

Ea 



1405.7570 

Ea 

31 + $t 


1402.7029 

Eg 

>31 + $t 


CMS-PAS-SUS-13-002 

uni 


Table 1; LHC searches used in this paper to test the viability of the simplihed models. 


generated HepMC event hies to estimate the efficiencies, e, of the signal regions dehned in 
all the ATLAS and CMS analyses that will be used in this work. The application examples 
and validation of Atom is found in [m \TE[ [TB] . We have implemented the CMS dilepton 
analysis in Atom and validated it using the cut-how tables given by the CMS collaboration 
based on the bi —)■ bx 2 bi~^i~Xi simplihed model. The comparison in the number of 
expected signal events calculated by Atom and CMS is shown in Appendix We also 
cross-checked some of the analyses with another simulation tool CheckMATE im. 

From the obtained cross section and efficiency, the SUSY contribution to the CMS 
dilepton excess is calculated as Nu = Ueff ■ C, where the ehective cross section, cJeff, is 
dehned as the cross section after the event selection: cTefr = e ■ Uprod- For the other ATLAS 
and CMS analyses the 95 % CL upper limit on aes, o'ul, is reported for each signal region 
by the collaborations. We dehne a useful measure for exclusion hj R = cXefr/cruL- If 
i? > 1 is found for one of the signal regions, the model is likely to be excluded, although 
one needs to combine all the signal regions statistically to draw a dehnite conclusion. 
However, we do not attempt to combine these signal regions because there are non-trivial 
correlations among them which originate from the uncertainties on e.g. the jet energy 
scale, the lepton efficiency and luminosity, and it is not possible for us to combine the 
signal regions correctly. Instead, in the next section we will look at the exclusion measure 
R individually to understand which signal regions are sensitive to the model points. 

In Table we list the ATLAS and CMS analyses we consider in this work. We include 
the multi-jet [HI |T9] and di-6 jet I2D] analyses, jets plus single EH or two lepton HlEiEa] 
(including same-sign (SS) dilepton [23]) analyses [22] and multi-lepton analyses [26l [271 
[ini|28]. In the next section we investigate whether the SUSY models can £t the CMS 
dilepton excess taking the constraints from these analyses into account. 
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Figure 4: Signal rate and i?-values for the squark models. The left panel presents the 
intermediate slepton, the right panel the off-shell Z scenario. 


5 Results 

5.1 Squark scenarios 

In Fig. I^we show the results of our numerical calculation for the squark scenario. In the 
plots the black curves represent the SUSY contribution, Ngg, normalised by the best fit 
value 130. The green bands correspond to the la region of the £t. In the same plots we 
show also the exclusion measure, R, for a few signal regions that are particularly sensitive 
to the models. The region where any R is greater than 1 is strongly disfavoured. 

In the left panel of Fig. [^we show iV^^/130 and R as functions of rriq for the squark with 
on-shell slepton model. We £x m^o = 495 GeV and m^o = 416 GeV. For these masses, 
there are no constraints from the chargino-neutralino direct searches. The right-handed 
slepton mass is fixed at 450 GeV such that iriedge in Eq. Q is 78 GeV, which is the optimal 
value for the GMS dilepton excess. 

As can be seen, this model can fit the excess only in the region where rriq < 650 GeV. 
However, this region is strongly disfavoured by the LI 10 signal region (shown in the blue 
curve) in the ATLAS stop search [25]. This signal region requires the same hnal state 
(2j as the GMS dilepton analysis, in particular an OS lepton pair with px > 25 

GeV and at least two jets with px > 20 GeV. The condition mx 2 >110 GeV is also 
imposed, which is very effective to reduce the ti and WW -|- jets backgrounds. One can 
see that the sensitivity of this signal region decreases as the increases due to the 
reduction of the production cross section. Nevertheless, the signal rate in the dilepton 
excess also decreases in the same way since these analyses employ similar event selection. 
Gonsequently there is no region in the plot where the SUSY events can fit the dilepton 
excess avoiding the exclusion from the other searches. This conclusion is robust against 
our assumption on the branching ratios, Br{qx q + Xi^/x^/Xi) = 10/5/85%, since the 




























LI 10 signal region constrains the same channel as in the CMS dilepton analysis. 

One can also see that in the rriq > 680 GeV region, the 2jm signal region in the ATLAS 
mnlti-jet search (TH] becomes sensitive and rules out the model points. This signal region 
is characterised by the requirement of at least two jets with px > 130 and 60 GeV and 
a moderately large effective mass, mefr = IPri’l + t^T > 1200 GeV, where is the 
Ath high pt jet with px > 40 GeV. The events with an electron or muon with > 10 GeV 
are rejected in this analysis. The 2jm signal region targets the qq —)■ topologV 

which is indeed the dominant event topology in this model since Br{qxi —)■ qx\) = 100% 
and Br{qx —)■ qxi) = 85 %j^ Due to the harsh cut on the rries, the 2jm signal region is 
sensitive to the models with large mass gaps between q and This is the reason why 
the sensitivity increases as rriq increases until the point {rriq ~ 850 GeV) at which a rapid 
degradation of the squark production cross section hnally turns the sensitivity down. 

In the right panel of Fig. we plot the V«/130 and R as functions of rriq for the squark 
with off-shell Z model, where we hx m^o = 478 GeV and m^o = 400 GeV so that rUedge 
in Eq. Q is 78 GeV. One can see that the SUSY contribution is too small to account for 
the dilepton excess, whilst this region is severely constrained by the 2jm and 4jl signal 
regions in the ATLAS multi-jet search [19]. Gompared to the on-shell slepton model, the 
rate of an OSSF lepton from a squark cascade decay is small: Br{qx —)■ qx^) ■ Br{x 2 
Z*x\) ■ Br{Z* —)■ £+£“) ~ 0.33 ■ 1 ■ 0.06 ~ 2 %, though we took a maximal value 33% 
for Br{qL —)• gy®) assuming xl Xi fo be Wino-like. Instead, xl Xi have large 
branching ratio to hadronic modes via Z* and W* which makes the off-shell Z model more 
prone to be excluded by the ATLAS multi-jet search [T9| compared to the on-shell slepton 
model due to the lepton veto cut in the analysis. The 2jm signal region constrains mostly 
QrQr QXiQXi topology and the sensitivity peaks around rriq ~ 900 GeV with m^o = 400 
GeV, similarly to the on-shell slepton model. On the other hand, the 4jm signal region 
requires at least 4 jets {px > 130, 60, 60, 60 GeV) and looks at the jets not only from the 
squark decay, q —)■ jx {x = Xi^X 2 but also from hadronic xf ^md X 2 decays and 

initial state radiation. Due to the milder cut meg > 1000 GeV, the sensitivity peaks at a 
much lower squark mass. 

We conclude that for the squark models it is very difficult to ht the observed GMS 
dilepton excess if the ATLAS stop search ESI and the ATLAS multi-jet search [T9| are 
both considered. 

5.2 Sbottom scenarios 

In this section we present the results for the sbottom scenarios, starting with the on-shell 
left-handed slepton model. In Fig. we show Na/lSO and R of the most constraining 

® We note that in m ATLAS does not exclude the region where m^o = 416 GeV in the squark- 
neutralino simplified model. We, on the other hand, exclude this neutralino mass for a certain range of 
the squark mass (see Fig. (left)). This is because our squark production cross section is larger than 
the atlas’s value since we set the gluino mass at 3.5 TeV in which the contribution from the gluino 
exchange diagram is still sizeable. 
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milGeV\ mi[GeV\ 

Figure 5: Signal rate and i?-values for the on-shell left-handed slepton mediated sbottom 
models. 

signal regions as functions of . As discussed previously, we assume m^o < — rrit to 

avoid tops in the decay chains that would lead to more stringent constraints. Within this 
condition we examine four different mass gaps: Am = — m^o = 50, 90, 130 and 170 

GeV. The left-handed slepton mass is fixed at = m^ — 40 GeV and m^o is set for 
each combination of m^o and such that rUedge in Eq. (fl| is 78 GeV. The intermediate 
slepton can either be a sneutrino or a charged slepton and the branching ratio of X 2 into 
these two states is assumed to be equal. Therefore, only half of the produced y® decay 
into an OSSF dilepton and a x?- 

In Fig. I^we see that a good £t can be obtained for sbottom masses between 420 and 
520 GeV, depending on Am. However, these model points are strongly disfavoured by 
the LlOO and LllO signal regions of the ATLAS stop search [25]. The event selection in 
the LlOO signal region is very similar to the LllO signal region which we briefly described 
in the previous subsection. The difference is that in the LlOO signal region the lepton p^ 
requirement is raised to (pt^Pt) > (100) 50) GeV and mT 2 > 100 GeV is imposed. As the 
lepton pt requirement is raised with respect to LllO, LlOO is especially sensitive to larger 
mass gaps. 
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Figure 6: Signal rate and i?-values for the on-shell right-handed slepton mediated sbottom 
models. 

In Fig. I^we show the contribution to the excess and the constraints from other searches 
in the on-shell right-handed slepton model for the four different Am, similarly to Fig.[^ In 
this scenario x° decays into an OSSF dilepton and a Xi with the branching ratio of 100 %. 
As can be seen, the results are similar to the left-handed slepton case and the region where 
the model gives a good fit is strongly disfavoured by the LllO and LlOO signal regions of 
the ATLAS stop search. The similarity of the results amongst the left- and right-handed 
slepton scenarios can be understood because LllO and LlOO constrain the same hnal state 
(2j -\- 2£ -f- as that is targeted in the CMS dilepton analysis and the kinematics of 
the dilepton events are similar between these scenarios. We conclude that it is difficult 
to attribute the CMS dilepton excess to the sbottom with on-shell slepton models if the 
constraint from the ATLAS stop search iza is taken into account. 

We now turn to the sbottom with off-shell Z models. The first model we investigate 
is the right-handed sbottom model where X 2 Xg are mostly Higgsino-like and Xi is 
mostly Bino-like. In this model the masses of three lightest neutralinos are calculated from 
the parameters, /i. Mi and tan/3, fixing M 2 at 3.5 TeV. Since we assume /x > Mi, we have 
m^o ~ ~ /i and m^o ~ Mi and both X 2 ^^nd Xg can contribute to the excess through 
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Figure 7: Signal rate and i7-values for the off-shell Z-mediated right-handed sbottom- 
higgsino models. 


their decays into an off-shell Z boson and a Xv The decay rate of the sbottom into the 
Higgsino states is dictated by the sbottom-bottom-Higgsino coupling which is proportional 
to tan/3. In order to have a large signal rate, we take tan/3 = 50 in our numerical scan. 
We again examine four different mass gaps Am = 50, 90, 130, and 170 GeV between the 
sbottom and X 2 - To this end we vary fi such that X 2 takes the desired mass set by Am. Mi 
is chosen such that m^o — m^o = 70 GeV. The mass of the lightest sbottom is calculated 
from given parameters hxing the left-handed third generation squarks mass, at 1.5 
TeV. A table with paramter values for each model point can be found in the Appendix 
In Fig. [^we show our results again in terms of A^^/130 and R. First we note that the 
strong constraint from LlOO and LllO observed in the on-shell slepton models is relaxed. 
To understand this we compare the distributions of mrp 2 , a kinematical variable used both 
in the LlOO and LllO signal regions, between the on-shell right-handed slepton (blue) and 
off-shell Z models (red) in Fig.at similar mass spectra. We take (m^^,m^o) = (400,230) 
GeV and hx m^o such that ruedge — 78 GeV for both models. For the on-shell slepton 
model we take m^ = 190 GeV. In Fig. [^we see that the off-shell Z model tends to give 
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Figure 8: Histogram of m'r 2 -distribution for sbottom production for both the off-shell 
Z and on-shell RH-slepton mediated case. The sbottom mass is 400 GeV and m^o = 
230 GeV. For the off-shell slepton mediated case we have = 190 GeV and m^o = 
151 GeV. The black vertical line indicates the cut for the limiting signal region SRA 
mGT150. 

smaller mT 2 compared to the on-shell slepton model. The solid (dashed) vertical black 
line represents the event selection cut on the mT 2 variable employed in the LlOO (LllO) 
signal region. As can be seen, the off-shell Z model is less sensitive to the the LlOO and 
LllO signal regions than the on-shell slepton model. 

What can also be seen from Fig. [^is that for Am > 90 GeV the SRA mGT150 signal 
region in the ATLAS di-bottom analysis 1201 is constraining and most of the preferred 
region of the dilepton excess is indeed disfavoured by this signal region. This signal region 
looks for two energetic 6-jets with pt > 130 and 50 GeV in events with tuct > 150 GeV[^ 
and $rp > 150. Events containing an electron [p^ > 7 GeV) or a muon {pj- > 6 GeV) are 
rejected in this analysis. This signal region is more constraining for larger Am because 
the event selection requires two energetic 6-jets. 

For Am = 50 and 90 GeV we hnd the regions where the observed excess can be 
explained at 1-a level without R > 1 from other searches. This result is consistent with 
the Endings reported in [6]. However these regions are already in tension with other 
searches. In particular the ATLAS stop search [23] and the ATLAS di-bottom search [23] 
give i? < 1 in these regions. 

In Fig. lolwe show A^^/130 and R as functions of m^^ in the left-handed sbottom model 
where X 2 iXi") is assumed to be Wino-like and decays predominantly to an ofi-shell Z 
(W) and a Bino-Iike x?. We again show the results for four difierent mass gaps and 

^ mcT = \j(^T “t -E't ~ (Pt ~ Pt )^> where Et and are the transverse energy and the transverse 
momentum vector, respectively. 
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Figure 9: Signal rate and -R-values for the off-shell Z-mediated left-handed sbottom model. 
Pure sbottom production is indicated by dashed lines and combined sbottom and stop 
production by solid lines. 


hx m^o = m^o — 70 GeV to fit the central value of the counting experiment. As we have 
mentioned in sectionjsl we assume the presence of the top squark, ii, with mi^ = The 
solid curves in Fig. [9 represent the results with both hih\ and iii\ production processes. 
To see the impact of the tilde tit\ production on the result, we also plot the contribution 
to and R from hih\ by dashed curves. 

One can see from Fig.j^that for Am = 50 and 90 GeV the model is strongly constrained 
by the SL5j signal region in the ATLAS jets plus 1-2 lepton analysis m This signal 
region requires a soft single electron (muon) with px € [10,25] ([6,25]) GeV and veto 
additional electron (muon) with px > 10 (6) GeV. It also requires > 5 jets with px > 
(180,25,25,25,25) GeV. The SL5j signal region is more sensitive to the iit\ topology 
where one of the stops decays hadronically G —)■ foyf —)■ hW*^ —)■ bqq'x'l and the other 
decays leptonically ii —)■ bxf —)■ bW*x^ —)■ biux^, because event selection requires a single 
lepton. We also note that the SL5j signal region becomes less sensitive for larger Am 
because the leptons from the stop cascade decay chain are boosted in this case and do 
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Figure 10; Variation of bi branching ratio into X 2 off-shell Z mediated left-handed 
sbottom scenario with no stop production. The color indicates Nu/130 and the black 
curves are lines of constant i?max- Only large wino branching ratios can provide a good fit 
to the excess. 

not pass the low pt requirement (< 25 GeV) efficiently. However, for larger Am the SRA 
mCT150 signal region becomes constraining. In particular the preferred region of the 
dilepton excess is disfavoured by this signal region at Am = 170 GeV. 

As a result we hnd a good £t to the dilepton excess at Am = 130 GeV and m^^ G 
[350,400] GeV, although this region is already in tension with the SL5j signal region in the 
ATLAS jets plus 1-2 lepton analysis. In addition, let us note that there is an additional 
constrain on the Xi ~ G mass plane from GMS single-lepton analysis [29], which is not 
included in our analysis. This analysis does not use the cut-and-count method but rather 
uses a BDT multivariate method, which prevents us from implementing this analysis. 
While recasting this analysis is out of the scope of this work, it is worthwhile to deduce its 
constraint on our models. Specihcally, the exclusion contour on the Xi~ti mass plane with 
chargino mass hxed at m-± = 0.25 rrti^ -1-0.75 m^o in the GMS analysis is most relevant to 
the allowed parameter space in our study ~ 0.3 m^^ -|- 0.7 m^o). At m^^ ~ 380 GeV 
the GMS analysis excludes m^o < 200 GeV, whilts m^o = 180 GeV at m^^ = 380 GeV in 
the bottom left plot (Am = 130 GeV) in Fig|^ 

It is interesting to note that the difference between the black solid and black dashed 
curves are small, whereas the difference is large amongst the blue solid and blue dashed 
curves in the bottom left plot (Am = 130 GeV) in Fig This means that the bibl pro- 
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duction gives the main contribution to the dilepton excess, while the model is disfavoured 


the contribution to the dilepton excess and the constraint from other searches in the 
versus BR{bi —)■ 6x2) plane concerning only the bibl production. In this study we assume 
BRipi —)■ 6x1) = 1 ~ BR{bi —)■ 6x2)- The region is divided into 3 colours, red, green 
and blue, which correspond to under, good and over £t of the dilepton excess, respec¬ 
tively. The R value of the most constraining signal region is shown in the black contours. 
One can see that a good fit is found for G [340, 380] GeV and BR(b —)■ 6x2) ~ 0-8 
without having R > 1 from other searches. Within our exploration we did not find the 
models where the sbottom is mostly right-handed and BR{b —?■ 6x2) ~ However, this 
result indicates that models that have a large cross section of the topology equivalent to 
bi —)■ 6x2 bZ*Xi can in principle explain the CMS dilepton excess avoiding constraints 

from other ATLAS and CMS direct SUSY searches. 


mainly by the additional titl production. Before concluding our study we show in Fig. 10 


6 Conclusions 


One straightforward supersymmetric interpretation of the observed dilepton excess by 
CMS lU |2] is the cascade decays of light-flavour and bottom squarks. In this paper, we 
studied and tested the viability of promising SUSY models by deriving constraints on these 
from various direct SUSY searches using the automated simulation tool Atom. 

In order to obtain a contribution to the dilepton excess from SUSY events, we con¬ 
sidered the decay of the second lightest neutralino, X2) either an off-shell Z or an 
intermediate on-shell slepton. The X2 itself arises from a light-flavour squark or sbottom 
decay. We investigated in total six possible simplified models, see hgures and 

We found that all of these models are already in strong tension with the experimental 
data once we demand a good £t to the dilepton excess. In particular strong limits arise 
from an earlier neglected ATLAS stop search [25] with identical hnal state topology. This 
analysis alone rules out the interpretation of the excess in terms of an intermediate (left- 
or right-handed) on-shell slepton for both light squark and sbottom production, see left 
panel of Fig. Fig. and Fig. respectively. We showed that if multijet plus missing 
energy searches are taken into account, the off-shell Z scenario with squark production is 
strongly disfavoured and noted that it is not able to give a sizeable contribution to the 
dilepton signal region, as can be seen in the right panel of Fig. 

We conhrmed the result reported in ^ and showed that the right-handed sbottom 
model with Higgsino-like X2 X3 decaying predominantly into an off-shell Z can explain 
the dilepton excess at 1-a level, although the model is already in tension with the ATLAS 
di-bottom search and the ATLAS stop search, as can be seen in Fig. This tension 
can be ameliorated if the left-handed sbottom model with Wino-like X2 H considered. 
However, the left-handed stop is necessarily introduced in this model and that creates 
another tension with the ATLAS jets plus 1-2 lepton analysis as be seen in Fig. 

We also showed in Fig. that in a simplihed model that only has sbottom production 
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the dilepton excess can be explained avoiding constraints from other searches in the region 
where — m^o 130 GeV, mi ~ 360 GeV and BR{hi —)■ 6 x 2 ) ~ 0.8, although we did 
not hnd a corresponding model point in the context of the MSSM within our exploration. 

This results may indicate that a more non-trivial SUSY scenario should be considered to 
explain the GMS dilepton excess. 

Note added: Shortly after this paper was submitted to arXiv, GMS updated their 
result [ 2 ] and reported most of the excessive events are observed associated with at least 
one 6 -jet. This new information further disfavours the squark scenario, which does not 
change our conclusion. Shortly after the GMS update, ATLAS released their new analysis 
of the jets plus SFOS dilepton channel inni. They explicitly looked at the signal region 
employed in the GMS dilepton analysis and did not hnd any signihcant excess. This 
casts a doubt that observed dilepton excess is merely due to the statistical huctuation or 
background mismodeling. The next run of the LHG will provide a dehnitive answer to 
this question. 

A Validation 

Here we show the validation results of our implementation of GMS-PAS-SUS-12-019/1502.06031 [H 
[ 2 ] and the ATLAS stop search with two lepton hnal state [25] . 

The benchmark point considered in the GMS analysis has a sbottom of mass 400 GeV 
decaying via b —)■ with 100%. The second lightest neutralino then undergoes an off- 

shell Z decay with SM branching ratios. We show the good agreement between the GMS 
results and our implemented analysis in Atom in table There, we give the event numbers 
in the central and forward signal regions as quoted by the GMS collaboration and their 
ratio to our results obtained with Atom. 

Additionally, we provide validation results for the stop search because of the strong 
constraints that we derive from this analysis. The ATLAS benchmark scenario consists 
of a stop decaying to yj*" -|- 6 with 100% probability followed by a decay of xf via a W 
into Xi and Standard Model particles. We show our validation in table In this table we 
present event numbers for the same-flavour (SF) and different-flavour (DF) case as given 
by ATLAS and their ratio to our results in the column Atom/Exp. 


{mi,m^o) = (400,150) GeV 

Gentral 

Atom/Exp 

Forward 

Atom/Exp 

Ajets Y 2 (no $rp requirement) 

242.7 ±2.8 

1.04 

34.2 ± 1.1 

0.77 

Ajets > 3(no requirement) 

186.2 ±2.5 

1.09 

25.6 ±0.9 

0.76 

> 100 GeV(no Ajets requirement) 

152.5 ±2.1 

1.03 

19.8 ±0.8 

0.98 

> 150 GeV(no Ajets requirement) 

85.0 ±1.5 

0.93 

10.4 ±0.5 

0.87 

Signal region 

132.4 ±2.0 

1.031 

17.0 ±0.7 

0.937 


Table 2: Validation table for our implementation of the GMS-PAS-SUS-12- 

019/1502.06031 analysis [H [2] in Atom. 
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(mf,m-±,m^o) = (400,390,195) GeV 

SF 

Atom/Exp 

DF 

Atom/Exp 

A0 > 1 

1834.9 

1.09 

2390.1 

1.06 


1402.8 

1.07 

1800.5 

1.07 

mT 2 > 90 GeV 

396.5 

1.02 

500.0 

1.09 

mT 2 > 120 GeV 

211.8 

1.01 

284.4 

1.1 

mT 2 > 100 GeV, PTjet > 100 GeV 

21.7 

1.4 

35.0 

0.99 

mT 2 > 110 GeV, prjet > 20 GeV 

86.0 

0.95 

116.1 

0.89 


Table 3: Validation table for our implementation of the ATLAS stop search with two 
leptons [22] in Atom. 

B Parameter Values for pMSSM scan 

In table 1^ we give additinal pMSSM input parameters as well as the sum of the branching 
ratio of hi to X 2 , Xs- These points were used to scan the right-handed sbottom-Higgsino 
model. Calculation of the physical SUSY masses and branching ratios was done using 
SPheno [E] E] . 



Am{b,X2) =50 

Am{b,X2) =90 

Am{b,X2) =130 

> 

too 

0 

mi =280 

(229,157,0.48) 

(183,157,0.63) 

(140,157,0.69) 

(96,157,0.63) 

mi =320 

(270,215,0.48) 

(228,215,0.63) 

(183,215,0.69) 

(140,215,0.67) 

mi =360 

(307,263,0.46) 

(269,263,0.63) 

(228,263,0.69) 

(183,263,0.7) 

mi =400 

(345,310,0.46) 

(307,310,0.63) 

(269,310,0.7) 

(228,310,0.71) 


Table 4: Additional information about the right-handed sbottom-higgsino model. We 
give values for (p, bNi)) for each model point. All masses and mass 

differences are given in GeV. See text for more details. 
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